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' The Gutzwiller variational wave function is shown to correspond to a particular disentanglement 

, of the thermal evolution operator, and to be physically consistent only in the temperature range 

U <^ kT <C Ef, the Fermi energy of the non-interacting system. The correspondence is established 
' without using the Gutzwiller approximation. It provides a systematic procedure for extending the 

ansatz to the strong-coupling regime. This is carried out to infinite order in a dominant class 
of commutators. The calculation shows that the classical idea of suppressing double occupation is 
' replaced at low temperatures by a quantum RVB-like condition, which involves phases at neighboring 

sites. Low-energy phenomenologies are discussed in the light of this result. 
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I. INTRODUCTION 



, Variational wave functions are a highly select class of results in the physics literature. There are only five which 

are widely used: Hartree-Fock and BCS for weak perturbations of the Fermi sea Feynman's ansatz for the ground 
state of liquid "^He 01, Laughlin's wave function for the fractional quantum Hall effect 0, and Gutzwiller's ansatz 
' for the ground state of the Hubbard model |^ . The last two fall into the class of Jastrow wave functions "S] , one of 
^ which was also used to describe "^He and ^He-^He mixtures and which are currently considered generic for strongly 
Q . correlated problems. 

^ ' Of these, Gutzwiller's is by far the least understood. Its underlying physical idea is that electrons of one spin see 
those of the other as a 'smeared background' ^^J- This very drastic assumption is still not sufficient to provide an 
operational prescription, but is supplemented by another, the 'Gutzwiller approximation:' electrons of one spin see 
^ ] the others 'as if occupying a band of width zero' i.e. their mass is taken to be infinite. This prescription was 
•■ never given a justification from first principles in any finite dimension. It is true by construction when the number 
of dimensions approaches infinity 7], because the scaling of hopping overlaps, required to obtain finite results in that 
^-H , limit, makes all motion effectively diffusive. 

' The present work approaches Gutzwiller's wave function from a perspective not suggested by its variational origin. 
It turns out that it is based on a one-step Trotter decomposition of the thermal evolution operator, strictly valid only if 
the on-site repulsion U is much lower than the temperature. This insight provides a natural scheme for improvement. 
A direct implementation of it shows that Gutzwiller's prescription to remove double occupation is the first step in a 
■ transcedent series. When summed, it yields a new projector, which imposes a quantum condition with much stronger 
I selectivity than the one removing double occupation. In the physical subspace satisfying this condition, Gutzwiller's 
I . program may be carried over to the strong-coupling regime kT <^ t <^ U as well, where t is the hopping overlap. 

' It can also be shown that at least at the level of expectation values, the on-site interaction does not scatter out 
Ch of the new physical subspace. Unlike the requirement of no double occupancy, the quantum condition cannot be 
O factorized into commuting local terms, indicating that relative phases on neighboring sites play an important role in 
the realization of the insulating ground state. The arguments are limited to the immediate vicinity of half-filling, 
where the configuration space for processes not considered here is small. 



■ II. THE GUTZWILLER ANSATZ 

Take the Hubbard Hamiltonian on a square lattice, H = K + V , where K is the kinetic term and V = U rii-^nii. 
Define the operator /C by a factorization of the imaginary-time evolution operator, 

^-f3(K+V) ^ ^-i3V/2^-l3K^-l3V/2_ ^-^^ 

The main result of this section is that the Gutzwiller ansatz neglects the entanglement of K and V. To prove this, 
take JC — K and calculate the expectation with respect to any operator O: 

tr e>e-^^/^e-^^e-^^/^ _ e-'^^/^0e-^^/^e-^^ \P) 

tr e-l^^ ~ iTpP^ 
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Here |P) are momentum eigenstates, and \R) , \R') position eigenstates. Use has been made of the fact that K is 
diagonal in momentum, and V in position: V \R) — UDj^ \R), where Dn is the number of doubly occupied sites in 
configuration \R). Now perform the same calculation for the expectation value {g\ O \g) in Gutzwiller's wave function 
\g) = P{g) 1^), with \^) a non- interacting ground state and P{g) the Gutzwiller projector: 

^(.9) = n [1 - (1 - 9)n^^n,^] = \{ [e, + 1,t + M + .9^] , (3) 

i i 

where the hatted operators are projectors onto empty sites, sites occupied by a single spin (up or down), and doubly 
occupied sites, respectively Then 

{g\0\g) ^ 5](*|i?)(i?|P(5)OP(5)|i?')(i?'l*) 

[g^-+^-' {R\0\R')]{R'\^) : (4) 
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remembering that P{g) \R) — g^" \R)- Now observe that the non-interacting ground state \^) is itself a momentum 
eigenstate. The same term P = will dominate the sum over P in Eq. if one takes the temperature low enough. 
The expectation value Q then reads 

RR' 

where [3* is a particular value of the temperature, for which the non-interacting system is in its ground state, to any 
desired accuracy. 

We are led to the astonishing conclusion, that the result Q of this procedure is the same as taking expectation 
values with respect to Gutzwiller's variational wave function, Eq. The correspondence 

{g\ O \g) ^ tr Oe-^'^'^e-^' e-^'^/^ (6) 

between the two expressions (0J and Q is established simply by replacing 

g^e-P'^'\ (7) 

The fact that the denominator in Eq. ||2Jl was tr e^^^ and not (more logically) tr e~^-^e~^^ is counterpart to the 
fact that {g\g) is not normalized. Taking O = 1, it immediately follows that {g\g) ^ tr e~'^'^e~'^^/tr e~P^ under the 
above correspondence. 

The original calculation of Gutzwiller also contains a prescription to fix /?*, or g. This is to take O = H, the 
Hamiltonian, and obtain g variationally. However, it is difficult to imagine such a procedure to compensate for the 
steps which were taken to arrive at Gutzwiller's form, Eq. ||SJ). Neglecting entanglement to get from Eq. (1) to Eq. (2) 
requires U <C kT, or more precisely, Ut <^C (kT)^. This is the step usually made in the Trotter formula, for a single 
short 'slice' of the evolution integral, which is eventually taken to zero {i.e. the temperature to infinity). To single 
out the ground-state term in Eq. (2) and so obtain Eq. (jsj, requires, on the other hand, the temperature to be low, 
kT Ep, the Fermi energy of the non- interacting ground state. The two are possibly consistent only in the range 

U-^kT-^Ep, (8) 

which is not the strong-coupling regime kT <^t <^U , for which the Gutzwiller approach was intended. 

On the other hand, in practice the expectation value is usually calculated in the so-called Gutzwiller approxima- 
tion, so it is desirable to understand its effect on the above derivation. It was shown in Ref. |8J that the approximation 
amounted to replacing the configurational overlaps by a constant, 

(i?'|P)(P|i?)^l, (9) 



3 



where TV is the number of terms in the sum over configurations R, R' . Inserting this in Eq. (0) yields 



RR' 



which is the same as obtained from Eq. Q with Gutzwiller's approximation, without invoking the hmit kT ^ Ep. 
The physical role of the Gutzwiller approximation is clear now. Instead of neglecting the excited states in Eq. (j^J, as 
done above by going to low temperature, it neglects the difference between the ground and excited states. When all 
the terms {R'\P) {P\R) are replaced by the constant l/M, of course their thermal average in Eq. ||2J) reduces to this 
single constant term. One gets again the same result as if the ground-state term alone had been taken into account. 

To summarize, calculating the expectation value of any operator with respect to Gutzwiller's wave function is 
exactly equivalent to the following three steps, when calculating its thermal expectation with respect to e"^^^^^^: 
first (a) neglect the entanglement of the potential and kinetic terms, so that one can replace 

^-P{K+V) ^ g-/3y/2g-/3Kg-/3V/2^ 

forgetting all commutators of V and K; then (b) take only the ground-state term from Eq. finally (c) use a 
variational procedure to fix the left-over temperature-dependent parameter e~^ ^/'^ = irrespective of consistency 
with the previous two steps. Needless to say, the first two steps are themselves hardly justifiable in the strong-coupling 
limit [/ > t. 



III. LOWEST-ORDER IMPROVEMENT 



In the previous section, it was shown that the expectation values calculated with the Gutzwiller ansatz can be 
obtained in a thermal formalism which neglects the entanglement of kinetic and potential energy terms. This is true 
irrespective of the use of the Gutzwiller approximation, which itself amounts to replacing the thermal average in 
Eq. |(2Jl with a normalization constant. The whole procedure retains only a classical attenuation of double occupation, 
such that Eq. is approximately valid with K, = K, i.e. all quantum dynamical correlations, coming from the 
commutators, are neglected. The main subject of the present work is to investigate the effect of the neglected 
commutators systematically. 

Including the commutators amounts to adding quantum correlations to Gutzwiller's wave function, which should 
be present in the strong-coupling low-temperature state, kT <^ t < U. Technically, this boils down to finding a better 
expression for the operator /C in Eq. The reason IC K is that the commutator [V, K] is not zero. Explicitly, 

V'^oK^tU'^Y. («,,_. - n,- ..)'= (aL V + ' (12) 

<i,3> 

where t is the hopping overlap, and the operation o is a commutator, 

V''oK= [V, V"-^ o K] (13) 

with V'^ o K = K . The vanishing of the commutator is obviously consistent with the original 'smeared background' 
interpretation, rii^-a — > {n-a-}- This points the way to an a posteriori justification of the Gutzwiller ansatz (though 
not of the Gutzwiller approximation). One can claim to work in a physical regime where it is sensible to replace 
the number operators by the average occupation of a site, as it should be in a doped metallic state, away from the 
metal-insulator transition. There one can hope that the high-temperature decomposition, Eq. pi|l . may in fact extend 
to low temperature. 

In the remainder of this section, the effect of including a single additional commutator will be studied. It will be 
shown below that the commutator [V, K] itself does not contribute to JC, because of the symmetry of the decompo- 
sition so the lowest non-zero contributions to same order in f3 are [V, [V,K]] and [K, [V,K]]. When U ^ t, the 
limit of interest here, the first is more important. Retaining only this one term, 

1^ ^ [V,[V,K]] 



6 V 2 



E 

<i,j> 
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(ala^^ + a],a,^) = K + Ki, (14) 
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where the numerical factors will also be justified later, and 



(15) 



is equal to one if the hop changes the number of doubly occupied sites, and zero otherwise. 

How can one use this result to improve Gutzwiller's ansatz? Note that the norm of Gutzwiller's wave function can 
be written 



^j. g-/3V/2g-/3ifg-/3V/2 



(16) 



tr e-f^K 

under the formal correspondence of the previous section. Obviously, one can interpret this as 

Pig) ^ e-^^/\ (17) 
To confirm the interpretation, recall the well-known alternative form 9] of writing the Gutzwiller projector (O, 



P{g) = cxp 



(18) 



where t] = — \ng. But this is just the right-hand side of (|17|l . under the correspondence Q. 

Now, the additional commutator in Eq. (|14|l amounts to replacing K hy K + Ki in the numerator of (|16|l . Since 
we have decided not to include any additional commutators, we may rearrange terms at will, and write 



g-/3\//2g-/3ifi/2g-/3/fg-/3i^i/2g-/3y/2 

tr e-f^^ 



(19) 



for the 'improved' right-hand side of (|16(l . It is obvious how to write an improved left-hand side now. There should 
be an additional projector, sensitive to configurations in which a hop would change the number of doubly occupied 
sites. Calling it Pi{gi), one may write 



Piigi) 



Explicitly, 



-Pi(ffi) = exp 



.9i 



En^T _.„{a^. a. -\-a\ a. 



where rji = — Inyi, and the new variational parameter corresponds to 

Hence the (systematically) improved variational wave function is 

|g,5i)=Pi(.9i)P(g)|*). 



(20) 



(21) 



(22) 



(23) 



The practical aspects of evaluating the projector 1)21(1 are beyond this article. One immediately obvious technical 
complication is, however, theoretically significant: the projector (|21|l cannot be written in the product form ||2J). The 
reason is formally that the individual terms in the sum in Eq. (|21|l do not commute. Physically, this means that even 
the first quantum correction to Gutzwiller's program for the one-band Hubbard model already involves the relative 
phases of fermions on neighboring sites. The importance of phases was already noticed in Ref. fl^, where it was 
shown that if a projector can be factorized into commuting local terms, then it cannot produce physical insulating 
behavior. In the following, it will be argued that while the use of projectors to remove unwanted parts of configuration 
space may be perfectly valid, such projectors cannot be used automatically to define effective Hamiltonians. 



IV. EXTENSION TO ALL ORDERS 



In the previous section, it was shown that the first quantum correction to Gutzwiller's projection involves bond 
phases. This is a setback for the local approach, so one is naturally led to ask, how important that correction is. 
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What is the range of vaHdity, in temperature, of the original ansatz? As luck would have it, this question admits 
of a sharp answer, because the correlation embodied in Ki (or Pi) above can easily be studied to all orders in jS. 
Unfortunately, the answer is rather disappointing: the low-temperature regime is approached exponentially fast. This 
statement should be moderated insofar as we really have in mind the state near half-filling. As already noted, any 
argument establishing a 'smeared background' automatically justifies Gutzwiller's ansatz at all temperatures. On the 
other hand, in the vicinity of the metal-insulator transition, one expects an electron to be scattered many times for 
each step it takes: there is not enough propagation in space to average out the local density. Then (|12|l contains the 
dominant processes affecting kinetic motion near half-filling, when U > t ^ kT . 

The resummation of all commutators in (|12() does not introduce new correlations, because 

(rii - n.j,^af^ = {n^-a - nj-af = n^j,-a (24) 

is independent of k. However, it will show how fast (in temperature) they become important, if one can calculate 
how these commutat;ors enter JC. The linear terms in a Baker-Campbell-Hausdorff (BCH) formula can be obtained 
by standard tricks ITl , to give 



sinh 



(3K ^ f3K 



K. = ' oK ^\\-^ coth ^ ) o y + Oiji^K^V'^) (25) 



for the case of Eq. (Some details are given in Appendix 1X1 ) The first term is now evaluated with the help of 
Eq. H12|) . remembering that the projector H24|) may be taken outside the Taylor series. One obtains 

/C = t ^ [1 + s{mn^,-o\ (ala^. + + 0(/3'i^') ^ + i^i + 0(/3^if 2). (26) 



where we have redefined K\. The point is now that the function 



sinh 

sm) = - 1 (27) 

2 

grows exponentially with pU . The interpretation H22|l of the variational parameter g\ should be replaced by 

g^^^-(it<mi-i_ (28) 

The projector P\ in Eq. becomes important at least as soon as 

s{m > 1, (29) 

when Ki and K become competitive. This is a much sharper condition than intuitively expected. On the other 
hand, it is also non-linear, so the correction in (|26|l is only 4% for kT = U, and the condition 129|l is first satisfied for 
kT « J7/4.4 . One may thus replace the lower limit of validity of Gutzwiller's ansatz by C//4 ^ kT (say), but that is 
obviously not essential. 

Comparing the interpretations |71l of Gutzwiller's parameter, and H28|) of gi, more can be said. When s{(3U) > U/t, 
the new projector Pi becomes more important than Gutzwiller's P. Thus the very-low-temperature regime kT U is 
completely dominated by Pi, unless t/U is exponentially small, which is not normally the case. (The limit C//t — > oo 
is discussed later on.) For kT <^ t < U , the wave-function Pi{gi) |^) represents reality much better (exponentially 
better, to coin a phrase) than P{g) |^). The strong-coupling limit is denoted t < U rather than t U to avoid 
confusion with the limit U/t oo, because we need s{(3U) U jt when kT ^ U . However, U jt is always taken 
to be sufficiently large to relegate the neglected commutators to weak perturbations. Note that, since s{f3U) rises 
exponentially, such a regime is easily achieved. For example, for kT = t/10 <t^t <U = lOt, one gets 

s{l3U) = s(lOO) - 10^^ > U/t = 10. 

The formal reason for the overwhelming dominance of Pi at low temperature is that the iterated commutator 
of Hubbard's contact interaction is non-zero to all orders. By comparison, for the harmonic oscillator already the 
third iterated commutator vanishes: [x^, ~ x^, and similarly when x and p are interchanged. If all possible 

commutators are arranged in a table, such that the {n,m) cell collects those of order V^K"^, this table is tridiagonal 
for the harmonic oscillator, while the expression H12|l gives the first column for the contact interaction. 
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The final question in this section is about the relevance of the other commutators in the above-mentioned table, 
of which so far only the first column was treated. At fc-th order in /3, commutators contribute which are of order 
j^myk-rn ^ ^rnjjk-m ^ rpj^g -jjj.g^ column, summcd Completely by the hyperbolic sine in Eq. H25|l . contains all terms 
with m = 1, i.e. like tU^~^ , while the other columns refer respectively to 1 < m < A;. Thus at any given order in /3, 
the term contained in the hyperbolic sine can be made to dominate those left out simply by increasing U jt. However, 
one should not hastily conclude that projectors generated by these terms are dominated by P\ at low temperature 
in the same sense as Gutzwiller's Pig). Like Pi, they are expected to have hyperbolic terms in (3U in the exponent, 
where P is only linear in PU. Thus they should be relatively mildly suppressed with respect to Pi, by a power of the 
ratio t/U. However, since suppression by Pi is the strongest, the configuration space determined by this projector 
alone is the largest one that needs to be taken into account, at least in the vicinity of half- filling. 

V. THE PHYSICAL SUBSPACE 

Looking back at the wave function (|23|) . in the light of the underlying dependence H28|l of gi with temperature, one 
may well wonder: what will survive such a projector at low temperature? Even if a wave function has a very small 
component to which Ki in (|26|l is sensitive, the hyperexponential suppression by Pi will annihilate it. The way out 
is obvious: the physically admissible space is the null-space of the operator Ki, Ki 1$} =0, or equivalently 

[V, [V, K]] 1$) - tU' + 1*^) = 0- (30) 

<i,j> 

This gives a precise meaning to the 'smeared background' condition at low temperature and near half-filling, where 
it cannot be trivially satisfied by the replacement n^.-o- — > 0. A more careful formulation of the same idea is 

exp [~(3{K + Ki)] 1$) = exp (-/3if ) |$) . (31) 

(The first of these conditions implies the second, but not vice versa.) However, Ki dominates K at low temperature 
because of the relative factor s{pu), so one expects the physics to be contained in Eq. H3U|) by itself. One may be 
tempted to add the condition 

V\<f)=0, (32) 

which is the no-double-occupancy constraint, but that is not warranted: for example, should the conditions H3U|) 
and (|32|1 turn out to be incompatible, the discussion in the previous section shows that the system will choose (|3Uf) 
at low temperature. By the same token, the subdominant correlations, coming from the neglected commutators, may 
have a say in which combination of the states |3>) turns out to be the ground state, but they cannot enlarge the 
physical subspace any more than the no-double occupancy constraint: Pi acts too stringently for that (if it does 
not, decrease the temperature, and/or increase the ratio U/t). If we decide to neglect all commutators containing at 
least two if' s, such as [K, [V, K]], because they are suppressed by at least t/U, the thermal expectation value in the 
physical subspace may be written at low temperature 

tr Oe-^'^/^e-^^+^'^e-^^/^ = ^ ($1 e-'^^/^Oe^'^^/^e-''^ [$) , (33) 

because the wave-functions which do not satisfy H30|l have been eliminated by the projector Pi, while the admissible 
ones allow the simplification H31() . This has formally the same structure as if we had made the high-temperature 
disentanglement Hll|) . i.e. used Gutzwiller's scheme. However, the underlying physical regime is at low temperature 
and strong coupling, ensured by the requirement (|30|l on the states |$). 

liU/t ^ oo, the dominance of Pi over Gutzwiller's P cannot be established quantitatively. This does not invalidate 
the former reasoning, but merely opens the possibility that the physical subspace is further reduced in the calcula- 
tion (|33|l . However, there is a qualitative argument that Hubbard's interaction leaves the physical subspace invariant. 
Namely, the identity 

VAV = i {V'^A + AV^ - [V, [V, A]]) (34) 

holds for any operators V and A. Now take V to be Hubbard's repulsion, and A = [F, [y, if]], c/. Eq. IpHljl . Let |$) 
be a state satisfying Eq. (|30|l . i.e. A\<^) = 0. Then the state V |$) also satisfies Eq. (|30|) . at the level of expectation 
values: 

{<^\VAV\^) = -i ($1 [V, [V,A]] 1$) cx ($1 = 0, (35) 
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because [V, [V,A]] oc A, by virtue of (|12() and (|24(l . In other words, Hubbard's repulsion does not on the average 
scatter out of the null-space defined by Eq. (|30|l . This is strong indication that the latter has been correctly identified 
as the physical subspace. 



VI. EFFECTIVE HAMILTONIANS 



In the previous section, the physical subspace was found to have two important properties. One, it is on the 
average invariant to V. Two, within it Gutzwiller's disentanglement ('smeared background') holds, so that hopping 
proceeds by the bare operator K, cf. Eq. (|3f |l . This should be contrasted with the no-double-occupancy subspace, 
which is strongly perturbed by K, requiring the introduction of projected hopping 

^ = ^ ^ (1 - n^^-a) al„aj„ (1 - n^- _„) + h.c, (36) 

<i,j> 

to keep within it. Both schemes take into account that an electron cannot hop locally onto anything but an empty 
site, in which case it tunnels by the full overlap t. In the first one, the burden of accounting for these dynamical 
correlations is taken by the construction (|30|l of the physical subspace, while in the no-double-occupancy scheme it is 
carried by the projected-hopping operator K. 

To compare the two clearly, note that using the projected-hopping operator in the no-double-occupancy subspace 
amounts to the calculation of thermal traces of the type 

tr Oe-^''^\-^''e-^^/' (37) 
at arbitrarily low temperature. Now, the exact expression for the trace is 

tr 0e-'3(^+^) ^ tr Oe-^^/^e-'^'^e-'^^/^, (38) 

by the definition of the operator /C. As shown in the previous sections, the relevant part of K. at low temperature is 
K + Ki, which in turn reduces either to Ki in the unphysical subspace (since Ki ^ K there), or to K in the physical 
subspace, defined to be the null-space of Ki. Even if this separation into physical and unphysical subspaces turned 
out to be wholly misguided, still the fact would remain, that JC has no component of the type K. The neglected 
commutators cannot change this, because (a) they can be independently suppressed, simply by increasing the ratio 
U/t, and (b) they contain K at least twice, as in [K, [V, K]], so they generate three-site and spin-exchange correlations, 
which are absent in K. 

In fact even more can be said: neither the classical constraint H32|) . nor the quantum constraint (|30|l can be used 
to define an effective Hamiltonian. The reason is that they are both negative statements, excluding some unwanted 
correlations. An effective Hamiltonian, on the other hand, acts to build up desirable correlations, not to reject 
undesirable ones. For example, the Hartree-Fock Hamiltonian has the respective Slater determinant as its ground 
state. By contrast, the null-space conditions (|30|l or (|32|l give no indication, which combination of states satisfying 
them is quasi-stationary with respect to the Hubbard Hamiltonian. Having an effective Hamiltonian is equivalent to 
knowing the approximate (saddle-point) equation of motion within the physical subspace, which is much more than 
knowing which states are not in that subspace. 

There is a simple and rather dramatic way to emphasize that null-space conditions are not equations of motion — 
or, equivalently, that projectors are not effective Hamiltonians. Within its regime of validity, one should be able to 
use an effective Hamiltonian H^jj just as if it were fundamental, i.e. forgetting that its elementary degrees of freedom 
have an internal structure. In particular, the same H^jj should be used in real and imaginary time. This reflects the 
requirement that the effective equations of motion admit the ensemble construction, i.e. are thermalized in the usual 
sense. Now, what is the real-time analogue of K + Ki, Eq. (|25|l ? It is found by taking /3 — ir in the function s{pU), 

sin (^) 

s{zrU) = - 1. (39) 

2 

Note the different dependence on U. To drive the point home, consider the limit tU — * cxd. Then s ~> ~l, and 



K 



Ki ^ i ^ [(1 - ni^^„)al^a^^{l - n^, _a) + (h.c. 

[?^^.-<TaJ^a^■^7^J■ _^ + (h.c.) (40) 



<i.i> 



t ^ 
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because the projector which appears may be rewritten: 

1 + s{iTU) n,,j_„ ^ 1 - = (1 - ni-a){'^ - rij-a) + n-i-an-j.-a- 

The first term in Eq. H40() is just the projected-hopping operator K , while the second makes the whole expression 
particle-hole symmetric. The appearance of K under these conditions is direct evidence that it is not an effective 
Hamiltonian for the one-band Hubbard model. If it were, it could be used in imaginary time as well, while it was 
shown above that the correct imaginary-time expression is Eq. 1)26(1 . not K. 

The formal way out of this real/imaginary time conundrum is again to work in the subspace l|3U|) . since the 
kinetic operator then reverts to the microscopic (bare) K on both axes. However, that is not the same as having an 
effective Hamiltonian. The prescription (|33|) correctly transcribes Gutzwiller's calculational framework to the regime 
kT <^t<U, but this is at best the first step in finding the normal modes, still far from revealing them. 

VII. DISCUSSION 

The standard operator approach to an effective theory without double occupation is to find a similarity transfor- 
mation |9l.[T^ 

SHS-^ = H,ff + ..., (41) 

such that Heff does not couple to the doubly occupied states to some given order, and truncate the remainder, 
represented above by dots. The present paper observes that such a procedure does not of itself guarantee a physical 
effective Hamiltonian. In particular, it is shown that avoiding double occupation is a semiclassical perception of 
the electron's behavior at strong coupling, which does not carry over to low temperature near half-filling. Imposing 
null-space conditions does not imply there exist thermodynamically stable normal modes which satisfy them. If the 
approximate normal modes are known, it may indeed be possible to connect effective Hamiltonians with low-lying 
states written in projector form, as in the cases of Hartree-Fock and BCS. If they are not, using the connection by 
formal analogy runs into trouble, such as not having the same effective Hamiltonian on the real and imaginary time 
axes. In other words, while an equation of motion may be written as a null-space condition, [H — i?) j^*) =0, not all 
null-space conditions are admissible constraints to some given equations of motion. Put more simply, fixing both the 
force and the effective constraint due to that force generally amounts to overspecifying the problem. To be certain 
the two are compatible is almost the same as knowing the solution. 

The projected-hopping operator K is a case in point, since it appears in the literature in two different contexts. 
On the one hand, it may be derived by formal arguments based on (|41|l . guided by the wish to avoid some parts of 
configuration space. On the other, the same K can be obtained by physical arguments from a more general three-band 
model, where it is claimed to describe the propagation of excitations against the background of a particular mode, the 
Zhang- Rice singlet this case, the use of if as a true effective Hamiltonian depends only on whether the correct 

hierarchy of background, excitation, and thermalization time scales is established. The point is that one cannot derive 
an effective Hamiltonian without some image of a physical mode which becomes quasi-stationary under the action of 
the original one. The same is true in the variational context. For example, in his description of ■^He-'*He mixtures Q, 
McMillan checked that the Jastrow wave function for the ''He background reduced in the long-wavelength limit to the 
density oscillations characterizing Feynman's formulation 0. 

There is nothing wrong, in principle, with using projector language to guess properties of the solution, i.e. to 
try and delimit the physical subspace. This is the essence of Gutzwiller's program, which the present work expands 
systematically. It is shown here that the original program is (at best) consistent in the physical regime U ^ kT ^ 
Ep, and that it may be carried over to the regime kT <^ t < U formally very simply, by working in the null- 
space of the 'most troublesome' commutator, Eq. I|30|l . This was identified with the physical subspace in the low- 
temperature, strong coupling limit, by a projection argument: states which are not in this null-space were found to 
be hyperexponentially suppressed at low temperature, much more strongly than states containing double occupation. 
The question may of course be raised, whether keeping the simplification of Gutzwiller's disentanglement at low 
temperature is worth the price of working with the condition H30|l . In effect, the preservation of disentanglement 
has replaced avoiding double occupation as a guiding principle for the construction of the physical subspace. This 
requirement is at least consistent between real and imaginary time, but it remains to be seen whether it is compatible as 
a constraint with the microscopic on-site repulsion. Physically, it amounts to the conjecture that low-lying excitations 
in the Hubbard model near half-filling can be mapped onto an effective semiclassical gas. It is encouraging for this 
point of view that the on-site repulsion does not on the average scatter out of the new physical subspace. 

The condition (|30|l is the first precise quantum formulation of Gutzwiller's 'smeared background' assumption. Its 
most interesting aspect is the role of phases. Indeed the expression indicates that the admissible physical states 
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should be coherent, as suggested by the resonating- valence-bond (RVB) arguments of Anderson 'T4'|. This should 
be contrasted with the viewpoint, based on the no-double-occupancy condition (|32|l . that the system is insulating 
because the electron locally has difficulty overcomin g th e repulsion U. The latter leads to an essentially diffusive view 
of the Mott state, which was shown already in Ref. 'l(f\ not to give physical insulating behavior, precisely because it 
is incoherent, i.e. insensitive to relative phases on neighboring sites. The fact that the interaction loses coherence, 
[K, V] — 0, as soon as nij,cr — 0, then implies that the ground-state of the one-band Hubbard model is likely to 
be antiferromagnetic at half-filling, even at arbitrarily large U/t. This conjecture is specific to the one-band model, 
since riij,^ appeares in the commutator only because both sites, connected by the hopping, are subject to the local 
repulsion. 

While this article was under review, new evidence appeared that the no-double-occupancy condition is not satisfied 
at low temperature in the one-band Hubbard model 15|. It turns out that the upper Hubbard band participates 
coherently in the low-energy density of states at half-filling, even for U/t as large as 12. 

To conclude, it has been shown that the Gutzwiller variational ansatz is physically consistent only in the regime 
U ^ kT <C Ef, inapplicable to the Mott transition. In the regime U > t ^ kT, the no-double-occupancy condition 
is replaced by a quantum condition (|30|l . sensitive to local phases, which defines the physical subspace. Within 
this subspace, Gutzwiller's disentanglement scheme is recovered in the low-temperature, strong coupling limit as 
well. However, as always, one cannot find an effective Hamiltonian without knowing the dominant slow modes. The 
quantum projection H3(J|I is hopefully a step forward in understanding their microscopic structure, but is not an 
equation of motion. 



Acknowledgments 

Conversations with S. Barisic and E. Tutis are gratefully acknowledged. Thanks are due to D. Svrtan for helping 
with Ref. ,11], and to P. W. Anderson for pointing out Ref. Q. This work was supported by the Croatian Government 
under Project 0119256. 



APPENDIX A: DERIVATION OF THE FIRST TERM IN EQUATION ^ 

Let A and B be two algebraic indeterminates. Define 

X = In (6-4+^6-4) , y = In (e-^e^+^) . (Al) 

It may be shown that the parts linear in B of these expressions are 

- 1 1 - 
X=^^—oB = Xb,Y=^ __oB = yB, (A2) 

where the notation =b means 'equal up to terms linear in B,' and the circle operation is a commutator, like in H13|l . 
The linear parts are called Xg and Yg, as noted. Now 

In (e-4e2-4+2Se-A) ^ (^r^X) (e^-e^-) Xb + Yb, (A3) 

where the second step is legal because neither X nor Y contain terms of order zero in B, and the third is trivial, since 
both Xb and Yb are linear in B. Adding Xb and Yg, one obtains 

In (e-4e2-4+2Se-A) e^_e^ ^ ^ ^^^^ 

This gives the first term in Eq. putting A = -pV/2 and B = -(3K/2. 

For completness, here is the derivation of ljA2|l from Ref. First, introduce a convenient notation for left- and 
right-multiplication by A, 

LB = AB, RB = BA, (A5) 

so that, for instance, the commutator with A is written (L — R)B. Then one can write A^BA"^ — L^R'^B, whence it 
is trivial to show 

e^Be-^ = e^e-^B = e^-^B = e^oB, (A6) 
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reverting to the circle notation. In the same vein, 



{A + B)" =B A" + ^ A'^BA''-^-'' = A" + ^ L''R"-^-''B, (A7) 
SO that the factorization formula 

n-l 

{L-R)Y^ L'^R"-^-'' = L" - i?" (A8) 

yields 

y4o(A + B)" eeb = (A9) 

from which follows the useful expression 

Aoe^+^ =^e^B-Be^, (AlO) 

upon summation over n. 

To get (|A2|) . first note that X =b Xb (no zeroth-order term), hence e-^ =b 1 + Xb- Then 

AoXB=n Ae^ - e""^ = Ae^+^e-^ - e^+^e^^A 
= (Ae^+s _ e^+s^) e--^ =b {e^B - Be^) e'^ = e^Be"^ - B = (e^ - l) o B, (All) 

where (lAlOjl and (|A6p were used in succesion. This means that in order to get Xb itself, we need one less commutator 
in each term on the right-hand side. In the circle notation, this is just 

e-^ - 1 

Xb = — ^ o B, (A12) 
which is the first expression in (|A2|) . The second is obtained in exactly parallel fashion. 
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